K* resonance effects on direct CP violation in B — > ttttK 
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Abstract 

B ± -> ir+ir-K ± and B° — > 7r + 7r K° decay channels are analyzed within the QCD factorization scheme 
where final state interactions before and after hadronization are included. The if* (892) and ^(1430) 
resonance effects are taken into account using the presently known ttK strange vector and scalar form 
factors. The weak decay amplitudes, which are calculated at leading power in AQco/rrib and at the next- 
to-leading order in the strong coupling constant, include the hard scattering and annihilation contributions. 
The end point divergences of these weak final state interactions are controlled by two complex parameters 
determined through a fit to the available effective mass and helicity angle distribution, CP asymmetry 
and K*(892) branching ratio data. The predicted iQ^lTSO) branching ratios and the calculated direct CP 
violation asymmetries are compared to the Belle and BABAR Collaboration data. 



PACS numbers: 13.25.Hw, 11.30.Er 



I. INTRODUCTION 



In the Standard Model, it is known that CP violation is mainly predicted in weak decays 
because of the weak phase of the Cabibbo-Kobayashi-Maskawa (CKM) matrix In decays of 
B mesons, through electroweak interaction, one can calculate the matter anti-matter asymmetry; 
VF-boson exchange and large beauty quark mass, allow a systematic perturbative calculation in 
the QCD factorization formalism (QCDF) 2] where the final state interactions are the main source 
of uncertainty. It is the combined occurrence of a weak and a strong phase differences that lead 
to the observation of the CP violating asymmetry between the B — > tt^tt^K and B — > tt^tt^K 
channels. 

Electroweak decays of resonant and non-resonant mesons made of a qq pair are well described 
in QCDF. In this framework, there is no direct three-body factorization scheme that efficiently 
describes a three-body decay, hence a quasi two-body state has first to be built up. Here, one 
of the two mesons assumed to be a K* resonance can decay via a strong decay mechanism to 
a (irK) state. In Ref. the authors attempt to reproduce the ttK effective mass and helicity 
angle distributions. In that calculation, the weak amplitude relies on effective QCD coefficients 
describing the leading order contribution as well as the vertex and penguin corrections at the order 
of h.Qcujjn^. The K* resonances decaying into ttK are then modeled by the scalar and vector form 
factors 4J that correspond to the strong final state interactions after hadronization. Additional 
phenomenological amplitudes, represented by four complex free parameters and added to the QCD 
penguin amplitude, are fitted to mainly reproduce the B — > K*(892)tt branching ratio and the CP 
asymmetry of the recent Belle and BABAR collaboration data. Furthermore they also predict the 
B — > i^*(1430)7r branching ratio. Altogether, one obtains a fair description of the data for these 
three-body B decays. 

In the present work, one explicitly includes the hard scattering and annihilation corrections at 
the order of Aqcd/^6- These weak final state interactions based on phenomenological assumptions 
are controlled by the endpoint divergences related to the asymptotic wave functions. This approach 
reduces thus the number of free parameters to only two complex ones. 

In Sec. m we derive the three-body decay amplitudes for the B — > ttttK processes within 
the QCDF framework introducing quasi two-body states. Sections IIIII and |TV] provide all the 
details for the weak decay amplitudes calculated at next-to-leading order in the strong coupling 
constant and in the perturbative expansion of the short distance interaction for B — > ttK*(892) 
and B — > 7r 1^(1430). Section [V] lists all the numerical parameters employed and in Section IVT1 a 
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discussion follows the presentation of the significant results on branching ratios and asymmetries. 
Finally, Section IVIII concludes with a summary of our work and some outlook. 



II. THREE-BODY DECAY AMPLITUDE 

To analyze the B — > ttttK decay amplitude, one first evaluates the matrix element 
(ttM2 I Heff I B ) within the factorization hypothesis, 

(TTM 2 \U eff \B) oc (M 2 \s lu (l- 75 )d\0)(Tr\uY(l-J 5 )b\B), (1) 

with M 2 being either the vector if* (892) or scalar Xq(1430) resonance, H e ff is the standard 
effective Hamiltonian for B decay (see Ref. [3|). The vector If* (892) and the scalar Kq(1430) 
resonances are assumed to be (ttK) quasi bound states in P and 5 waves, respectively. Thus, one 
writes Q], 

(mrK\H ef f\B) oc ((7rK) 5> p|s7 ! ,(l-75)d|0)(7r|n 7 "(l-7 5 )6|5), (2) 
where ((ttK)s,p \ s~/ u (l — 75)^ | 0) is expressed as 
(ir( P7T )K( PK )\s 7u (l- l5 )d\0} 



,„2 _ m 2 

(j>K ~ Ptt) 
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f^^) + !±-< qv f^^). (3) 



In Eq. (|3j), q 2 with q = px + Pw is the invariant ttK mass squared, mx and denote the kaon 
and pion masses, respectively. The vector ff K {q 2 ) and scalar fo K (q 2 ) form factors are describing 
the final state interaction after hadronization. From semileptonic decays like r — > Kttv t and 
K — > nlvi, one can extract informations on these Kir scalar and vector form factors [4]. Analyticity, 
unitarity, QCD counting asymptotic rules allow one to relate scalar and vector form factors to the 
^(1430) — > irK and if* (892) — > ttK scattering amplitudes in the elastic and inelastic domains. 



All the details can be found in Refs. 



. [3] and [i 



[4]. The full amplitude for each wave is given by 



A 3 (B -> ttttK ) = A(B ttM 2 ) x T(M 2 Kir). (4) 

For the if* (892), the vertex function r(K*(892) -» Kir) associated with the B ttK* ttttK 
decay is written as 

r(K*(892) -+Kn) = -?- gzEl-EzEI f**tf), (5) 

qjK* e K , (892) -p B 

where fx* is the K* decay constant and e^-.( 892 ) • Pb = (mB/^q) A 1 / 4 (m 2 B , q 2 , m 2 ) , ps and mg 
denoting the B four momentum and mass, respectively. In Eq. (0), A(x, y,z) = (x + y — z) 2 — Axy 



and the moduli of the ir^ momenta are, 

1 

Ptt+ = TT" 

2q 

and, 



K - m w y 



1 

2q 



m 2 B -(q + m n ) 2 m\ - (q - m n Y 



The vertex function r(^(1430) -> Kir) is 

r(if5(1430) -> Kir) 



Ik. 



1 m | - ml K 2 

—9 JO W J 5 



(6) 



(7) 



(8) 



where /^-* denotes the JTq decay constant. Following closely [2j], one derives the QCDF decay 
amplitudes where the short and long distance contributions are factorized in the approximation of 
a quasi two-body state, 7rif*(892) or 7rET *(1430). 
The amplitude B n~K*°(892) is 



q=u,c \, 



5 qu AO) + 40) + *0K0) 



- ^ioO) + r ^*0)«sO)) + AO) + A,ewO) 



>, (9) 



and the 5° -> vr+K*-(892) amplitude, 



q=u,c 



5 qu a\(n) + a\{n) + r** (//)ag(/i) 



1 



+ afoO) + r x OKO) + AO) - ^A,ewO) 



■, (10) 



where the coefficients a^O) and /3 n O) are given in Eqs. CLE]) and ([20]). The A^ s) are product of 

CKM matrix elements, the r^ 2 (fi) the chiral coefficients and [i is the scale. 
1/1 /iQn\ 



The 5- -> 7r- J PC*°(1430), amplitude reads, 



q=u,c 



<VAO) + alO) -r x °(n)al(n) 



o a H)0) - r X ° OKO) + AO) + A,EwO) 



(11) 



while the B° -)■ 7r+.K"o~(1430) amplitude is 



.4(3° a**?-) = E 4 S) v 4 ^* 



q=u,c 



Squ a?0) + 40) - OKO) 



+ a ?oO) - r ?° OKO) + AO) - t:A,ewO) 



(12) 



ttK* 
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The chiral coefficients, r^* (p) and r^* ] (p), will be given in Eqs. (|24p and (|40p . 



For the JT*(892) resonance, the pseudoscalar- vector factor A n K* in Eqs. ([9|) and (fTOjl reads, 

•W ' ~^ ■2 (l< l ..,,^.p B F^(q 2 )f K ., (13) 



~ l e ^*( 892;i 



with the Fermi constant Gf = 1.16 x 10~ 5 GeV~ 2 and where the weak transition form factor 
F^^ n (q 2 ) will be given in Sec. El For the 1^(1430) scalar resonance, the pseudoscalar-scalar 
factor A^k* in Eqs. (jlip and (|12p . is 



^ = i ^(m| - m2)F ^(g 2 )/^. 



(14) 



In Eqs. (|9|)- (|T2|) . the CKM matrix elements are, 



X[ s) = V ub V: s = AX 3 (p-irj) A, 



\i s) = V rh V r 



\2 

A\ z I 1 — 4- 



f. E 



(15) 



where following Ref. [5|] the Wolfenstein parameters are, A = 0.814, p = 0.1385, rj = 0.358 and 
A = 0.2257. 

Since one assumes the dominance of the if* (892) and -K"q(1430) resonances in the description 
of the irK channel, the full amplitude A%(B — > ttttK) is built up on the P and S waves so that the 
differential effective mass branching fraction is [3J, 



dB(B — )■ ttttK) 1 q \p n +\ \p n - 



A 3 (B -> tt(itK)s) 



1 

+ 3 



A 3 (B -»• 7r(vrK)p) 



dq T B 4(2vr) 3 m| 

where r# = 1/tb is the P-decay width. The usual CP violating asymmetry parameter is 

. ,„ _ B(B ^ ttttK) - B(B ^ ttttK) 

A C p(B -> ttttJC) = —3 e r= =4- 

v y B(B -> vrvriO + B(S -»• vfTfK) 

In Eqs. ([9|)- (fT2]) . the a^(/z), involving the Wilson coefficients C n (p), are 
C„±i(m)" 



(16) 



(17) 



C„(/i) + 



+ 



Cf 



AirN r 



N n (M 2 ) + J*(Af 2 ) 



a s (p)C n±1 (p)V n (M 2 ) + 



^ 2 a s (p/2) 



C n±1 (p/2)H n (irM 2 



(18) 



with n E {1, 10}, and the scale is = m&, m& being the b quark mass. In Eq. (|18p . the color 
number is iV c = 3 and Cf = (N 2 - 1)/2N C = 4/3. The upper (lower) signs in C n ±\(p) apply when 
n is odd (even) and 



N n (M 2 



0, n e {6,8}, and M 2 = K* (892), 

1, else. 



(19) 



The Wilson coefficients, C n (fi), computed in the Naive Dimension Regularization (NDR) 
scheme Q], are taken at the scale nib for the vertex, V n (M 2 ), and penguin P^{M 2 ), corrections 
which involve only the fo-quark, whereas the annihilation, fi n (jxM2) and hard scattering, H n (M2), 
contributions are evaluated at the scale nib/ 2 since they involve the spectator quark. The strong 
coupling constants are a s (nib) = 0.224 and a s (nib/2) = 0.286 [jj. 

The annihilation term, f3 n {u), is given by 

\a , m [^(M)]7rM 2 B *M* , orA 
\PM\*M2 = A ' ( 2 °) 

where the factor, B w m 2 , ls the product of Gf by the B, tt and M 2 decay constants, 



G F 

BirM 2 = =R fBfnfM 2 , (21) 

with the upper sign if M2 = K*{892) and the lower sign otherwise. In Eqs. ([9])-(|12p. the tree 
annihilation component (at /x = nib/ 2) is (the upper-scripts I and F denote initial and final 

states), 

M^M 2 = %C2{li)A[^M 2 ), (22) 
while the penguin annihilation terms (at fi = nib/ 2) are 
_ Cj 



A/ 2 jy2 



C 3 (n)A[(irM 2 ) + C 5 (n)[Al(nM 2 ) + Af(7rM 2 )) + N c C 6 {ji)Ag(<irM 2 



[ fe 3,Ew(/i)t M2 = ^2 



C 9 (/x)A{(7rM 2 ) + C 7 {ii)(a{{ttM 2 ) + Af (ttM 2 )) + N c C 8 (fi)A^ (ttM 2 



(23) 



where the amplitudes A 1 -^ (7rM 2 ) are given in Eqs. ([36]) for the P-wave and ([50]) for the 5-wave. 



III. QCDF CORRECTIONS FOR B -> ttK*(892) 

The pion in the final state ttK*(892) is created from the transition B — )■ tt while the K*(892) is 
created from the vacuum; this mechanism is due to the structure of the four-quark operators in the 
heavy quark effective theory as well as the conservation of the flavor quantum numbers. Following 
Ref. {21, we only give the QCD corrections that appear in A{B — > tt K* (892)). 

Since the coefficients in the Gegenbauer expansion of the light cone distribution amplitudes 
(LCDA) are known with large uncertainties [2|, one here limits oneself to leading terms in this 



expansion for the tt and K* (892). The leading twist- 2 distribution amplitude is <J?(x) = 6x (1 — x) 
and the twist-3 two particle distribution is (p(x) = 1 and ip(x) = 3(2x — 1) for both ir and K*(892). 
The chiral coefficient for the vector meson if* (892), given at the scale fi is defined as 



„K* r ,,\ _ ^\f<? fx* 



r? (M) 



(24) 



m b (n) fx* ' 

where /^« is the transverse decay constant and where one has introduced the running meson mass 



square, replacing rn^, 8 g 2 \ by m 2 iK = q 2 . For a pion, the chiral coefficient reads 



rut 



m b (fi)m u (^) 



(25) 



with the u-quark mass m u . 



A. Penguin contributions 

The penguin contributions, P„(-K"*(892)), with the values n = 4, 6, 8, 10, required in Eqs. 
and ([TO]) , are as follows, 



4, m b 2 

3 ~u" 3 ~ Ga '*( 892 H s <?) 



8 7Y}/fo 4 

2 In — + 2 - G ! Jf*(892)(0) - G , ^»( 892 )(l) 



+ (C 4 (/i) + C 6 (/x)J x 



+ <2»(m) 

(«F - 2) Gx*(892)(0) - G^.( 892 )(s c ) - G K *( 892 ) (1) 

-*Jf(892)(a:) = 3, 

and C|g (//) related to the Qsg chromomagnetic dipole operator. Furthermore, 



4n F . m b 

in — 

3 n 

(892) ix) \ , (26) 



with 



(27) 



P|(iT(892)) = -^M 



Cl(/x) Gk*(892)(s 9 ) + C 3 (/i) 



C/r*(892)(0) + ^.0(892) (1) 



+ c 4 (m) + c 6 (aO 



(Wi? - 2) G^*(892)(0) + Gk*(892)(Sc) + Gr*(892) (1) 



P|(iT (892)) 



a, 



9vriV r 



dfji) + N c C 2 (fJ.) )G K 

*(892) 



(28) 



(29) 
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where a e = 1/129 is the electromagnetic coupling constant. Finally, 



Pf (if*(892)) 



9vriV r 



Ci(/x) + N c C 2 (v) 



4 ! m 6 2 



-SCf^/x)^ 1 ^^.^)^!. (30) 

In these equations, /i = mft and the number of active flavors is np = 5. In Eq. (|30p . C^(/x) is 
related to the electromagnetic dipole operator. The gluon kernel contributions are 



Gk* (892) (Sq 



5 2l7T 

1 1 S<7 = 0, 

3 3 q 
85 /- 4tt 2 



= 1, 



5 2 32 o 2 I 

- - -ms c + — s c + 16s c - -VI - 4s c 

4 



-17T 



2 arctanh (Vl — 4s c ) 



1 + 2s c + 24s^ 

2 



?7T 



2 arctan h (Vl — 4s c ) 



(31) 



and 



Cx*(892)(Sg) 



1, Sg = 0, 



4vr k 



-35 + 4V37T + — -, s Q = 1, 



(32) 



-12 sf 



2 arctanh (Vl — 4s c 



ITT 



36 S r 



+12 Vl - 4s c s c 2 arctanh (Vl - 4s c ) 



17T 



~t~ 1; *->q — ^C* 



In Eqs. ([3T]) and ([32]) . s 9 is defined as (rag/nib) 2 so that s g = for g = it, <i, s 9 = 1 for q = b and 
s q = s c for q = c. 

B. Vertex contributions 



In the B — > irK*(892) transition, the electroweak vertex, V n (K*(892)), receives a s (/i) correc- 
tions to all a1i{y) in the amplitude A(B — > irK*(892)), 



V n (K*(892)) = { 



' 12 1nf^^ -3wr- y, n£ {1,4,10}, 



(33) 



9-6tvr, nG{6,8}. 
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C. Hard scattering contributions 



Evaluated at the scale [i = m&/2, the hard scattering correction can be written as, 



HjnMo. 



B-kM 2 
A n M 2 



H n {nM 2 ), 



(34) 



where for the vector resonance, M 2 = K*(892), ^4,rK*(892) an d ^7rX*(892) are defined by Eqs. ([13 
and (f2"Tj) . respectively. One has 



H n (nK* (892)) = < 



3^ [r^ji) X H + 3] , n €{1,4,10}, 
0, n£{6,8}, 



(35) 



where As = 0.3 GeV is a hadronic parameter of the order of Aqcd ©]• In Eq. ([35]) . r^(fi) is given 
by Eq. (|25p and X// represents the end point divergence related to the soft-gluon interaction with 
the spectator quark. Its expression will be given in Eq. (f5Uj) in Sec. [V] 



D. Annihilation contributions 



The annihilation amplitudes cannot be derived from the QCDF approach so that they are 
model-dependent involving also a divergence parameterized by Xa (Eq. I[56|) ). Based on Ref. Q], 
the expressions for A I j (-irK*(892)) and Af (vrK*(892)), for j = 1 and 3, are, 



A i 1 (ttK*(892)) w 6vra s (/i) 
^(7riT(892)) ps 6to s (^) 

A%(ttK*(892)) ss -6vra s (^) 
with = m&/2. 



3 ( A' i - I + y ) + rf (^a " 2X A ) 



-3rf (/i) (xi-2X A -^ + 4) +rj(/i) ( *1 - 2X A + y 



3rf *(/i) (2X A - 1)(2 - X A ) - r£(/x) (2*1 - X A ) 



(36) 



IV. QCDF CORRECTIONS FOR B -> 71-7^(1430) 

We now turn to the B — > 7tA"q(1430) transition for which the a s (fj) corrections are all included. 
Here again, only the first non- vanishing leading term in the LCDA of the A"q(1430) are retained: 



<1>k*(1430)OzO = 6x(l -x)\l + 3< S(143 °V) (2x 



(37) 



where S^ ^ 1430 ^// = mj) = 5.26, and B^"^ 1430 ' 1 (// = m&/2) = 0.39 are the first non- vanishing 
Gegenbauer moment (for neutral scalar) evaluated at two different mass scales. The asymptotic 
form of the LCDA for the pion is 



$ n (x) = 6x(l - x) 



(38) 



The twist-3 two particle distributions are 



<Afq (1430) 0*0 = 1 and tp^x) = 1. 



(39) 



Similarly to the B — > itK*(&Q2) decay channel, the B — > 7t.Kq(1430) decay amplitude is fac- 
torized out into a product of a transition form factor B — > ir times a A^q(1430) decay constant as 
shown in Eq. (2) of Ref. [2]. 



The -KTq(1430) chiral coefhcient is given by: 



„tf *(1430) 



00 



2q 2 



(40) 



m 6 (/z) (m s (/x) - m u (n)) ' 

where m s is the strange quark mass. In Eq. (|40p . as has been done for the if* (892) meson (see 



Eq. (|24p . one has introduced the running meson mass square for the -K"q(1430) replacing 



*(1430) 



by m 



nK 



A. Penguin contributions 



f. Q 



From Ref. [2j, one can obtain all the penguin corrections P,? (if q (1430)), (with n = 4,6,8,10) 
for the B to pseudoscalar-scalar transition. One has, 



P 4 9 (A- *(1430)) 



4vrA^ r . 



Ci(m) 



4 m 6 2 

3 m — + 3 - G E- *(1430) (, s gJ 



C 3 (/*) 



3 ln ~^ + 3 ~~ G '^o(i430)(0) - G , A -»( 1430 )(1) 



x Cx*(1430)(0) - G\fC*(1430)( s c) - Gx*(1430)( 



+ (GiO*) + C r 6 (M) 

1 dx 



^ln^-^-2) 
3 /i 



2C £Mjf ^^(1430)^) , (41) 



with 



— $K* o{ um)(x) = 3B 1 °(v) + 3. 



(42) 
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Moreover, 



^(^(1430)) 



Ci(M) 



4 m 6 2 
-In — + - - G 
6 ii 3 



Gk"*(1430)(0) _ G\ftr*(i430)(l) 



A*(1430H s gJ 



+ [C i (ji) + Ce{fi) 



8 m 6 4 

- in 1 — 

3 [i 3 



4n F m b 

in — 

3 /i 



( n f - 2 )Cii-*(i430)(0) - ^.(1430) (s c ) - Grt*(1430)(1) 



2CSf(M)S, (43) 



Pf(K *(1430)) 



a, 



9vriV r 



and 



Pf (if *(1430)) 



9irN r 



diii) + N c C 2 {n) 



In — + - - G^*( 1430 )(s (? ) 



3C?» , (44) 



C^n) + N c C 2 (fi) 



3 ln ~ + 3 - G ^ *(l430)(s g ) 



SC^^^Y^^^O)^)!. (45) 



Comparing Eq. (|26|) with Eq. (|41|) and Eq. (|30|) with Eq. (|45|) . one can see that the formal struc- 
tures of P|(^o( 1430 )) and Pi ( K o (1430)) in terms of C n (fi), of gluon kernel functions, G M2 {s q ) 
and of LCDA, &m 2 { x )> where M 2 is now ^^(1430) instead of K*(892), are identical to those of 
P%(K*(892)) and Pi (K* (892)) , respectively. The gluon kernel functions, entering in Eqs. (|4"Tj) - 



5]), are 



KX, 



5 2i?r (/i) 



3 + - + 



2 

4vr 2 



85 /- 

T - 6 ^ + - 



155 



G\A*(1430)( s <j) — < 



5 2, Sf°*(/u) 4 

lns c H 1 K ' + - 

3 3 2 3 



0, 



36\/3vr + 12tt 2 



8 + 95f°*(/i) s c + 2 S + GSSf ^) 



,A* 



s q = l, 



306#f 5 (/x) S 3 _ 1^/1 - 4s c 1 + 2s c + 6(4 + 27Bf° {u,))s 2 c - 324B* 5 (p)a\ 



A'* 



,A ( ? , 



2 arctanh (\/l — 4s c ) 



^7r 



+ 12s" 



l + 3B^(/i)-|(l + 9flf 5 



18Bf 6 (/x)s 2 2 arctanh - 4s c ) 



27T 



(46) 



and 
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Gjf*(1430)( s g) 



( 16 2vr 

I i, 

9 3' 

-32 2tt_ 

~ + 7s 

16 2 r 

— (1 — 3s c ) - - lns c + (1 - 4s c ) 3/2 [2arctanh (Vl - 4s c ) - ot] 
,9 3 L 



Sg = 0, 



So = 1, 



1 Sg — S c 



(47) 



B. Vertex contributions 



The relevant vertex corrections, V n (K g (1430)), for n G {1,4,6,8, 10} are the following, 



F n (tf *(1430)) = < 
with /i = m^. 



— M _3i^__ + _(ii_6i7r)sf S (^) > n£ {1,4,10}, 
-6, nG{6,8}, 



(48) 



C. Hard scattering contributions 

From gluon exchange between the scalar 1^(1430) and the spectator u quark one derives, at 
/i = mft/2, the hard scattering corrections. One writes it as in Eq. ()34|) with M2 = -fTo(1430) and 
J 4 7r ^*( 14 30) and B nK * (1430) defined by Eqs. (Q3J) and (|2"T]) . respectively. Here 



3ms 



F n (^*(1430)) = { 



3(Bf 0, (//) + l)-r;(M)X H (^(M)-l) , nG{1.4,10}, 



0, nG{6,8}. 

As for the If* (892) (Sec. IIII Cl) the endpoint divergence is modeled by Xh- 



(49) 
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D. Annihilation contributions 



The weak initial and final annihilation amplitudes, A^Tri^Q (1430)) and AJ(ttKq (1430)) (with 
j = 1 and 3) at (J, = rrib/2 are calculated starting from Ref. B for B -4 vrifg (1430): 

A{(itK* (U30)) « 2vra s (/i) f 9flf 5 (p) (3X A + 4 - tt 2 ) - r£(/i) r* S (») *lj , 

4(7T# *(1430)) « 67ra s (/i) j^MB^O*) [x\ - 4X A + 4+ + rf° ( M ) - 2V A + yj | 
Af (vr^(1430)) « 6vra s (/i)X A |^( Ai )sf *(^)(6X J 4 - 11) - r?'(/i)(2X A - 1) j, 

(50) 

with an endpoint divergence (Eq. (|56p). These amplitudes will be then implemented in the 
b n (nK£(UW)) given in Eqs. (JSJ) and ([23]). 

V. INPUT 

A. Numerical parameters 

In this Section, one summarizes all the values of the parameters required for performing numer- 
ical applications. From Ref. B, the meson masses in GeV are 

m B = 5.300, = 0.139, m K = 0.493, m K * = 0.892, m K * = 1.430, m B * = 5.320. (51) 

The running quark masses (at = 4.2 GeV) in GeV are 

m b = 4.2, m c = 1.3, m s = 0.070, m u<d = 0.003, (52) 

whereas at rrift/2, one has in GeV {?]], 

m b = 4.95, m c = 1.51, m s = 0.090, m u4 = 0.005. (53) 

The meson decay constants in MeV are 

f B = 180 ± 40 fl, f K * = 218 ± 4 Q, U = 130 ±0.2 Q , fit* = 175 ± 25 Q- (54) 

The scalar meson decay constant fx*, which appears in Eqs. ©, (HU) and PI]) , does not, in fact, 
enter in our calculation as it cancels out in A%{B — > irirK) [Eq. fl3J)], in f3 n {i:M2) [Eq. (120]) ] and in 
H n (irM 2 ) [Eq. flM}]. 
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TABLE I: Real and imaginary parts of the leading order (LO), vertex, penguin and hard-scattering contri- 
butions to the short distance amplitude, a^(/x), for P and S waves. The scale /i = mi, except for the hard 
scattering where fi = rrib/2. 



P-wave 




LO 


Vertex 


Penguin 


Hard scattering 


Total 


a i w 
ai{n) 


(1.018; 0) 


(0.028; 0.014) 


(0;0) 


(-0.246; 0.317) 


(0.800; 0.331) 




(-0.031; 0) 


(-0.002; -0.001) 


(0.003; -0.014) 
(-0.002; -0.005) 


(0.018; -0.023) 


(-0.012; -0.038) 
(-0.017; -0.029) 


G§(m) 


(0;0) 


(0.0006; -0.001) 


(-0.007; -0.0009) 
(0.001; 0.011) 


(0;0) 


(-0.006; -0.002) 
(0.002; 0.010) 


a u s (p) 


(0;0) 


(-0.6; 1.3) x 10~ 5 


(-4.7; 0) x 10" 5 
(-0.3; 6.4) x 10" 5 


(0;0) 


(-5.3; 1.3) x 10" 5 
(-0.9; 7.7) x 10" 5 


aio(M) 


(-0.0014; 0) 


(0.0014; 0.0007) 


(0.0002; -0.0001) 
(0.0002; -0.0001) 


(-0.013; 0.016) 


(-0.012; 0.017) 


>S-wave 


o¥(m) 

a\(n) 


(1.018; 0) 


(-0.016; 0.089) 


(0;0) 


(-0.151; 0.184) 


(0.851; 0.273) 




(-0.031; 0) 


(0.001; -0.007) 


(0.023; -0.017) 
(0.039; 0.036) 


(0.011; -0.014) 


(0.004; -0.037) 
(0.021; 0.016) 


g§(m) 


(-0.039; 0) 


(-0.0004; 0) 


(-0.003; -0.014) 
(-0.006; -0.004) 


(0;0) 


(-0.042; -0.014) 
(-0.045;-0.004) 


a u s (p) 
ag(/z) 


(44; 0) x 10~ 5 


(0.4; 0) x 10~ 5 


(4; -10) x 10~ 5 
(2; -5) x 10~ 5 


(0;0) 


(48; -10) x 10- 5 
(46; -5) x 10~ 5 


a"o(M) 
a io(M) 


(-0.0014; 0) 


(-0.0008; 0.005) 


(0.0015; -0.0001) 
(0.0016; 0.0002) 


(-0.008; 0.009) 


(-0.009; 0.014) 
(-0.009; 0.014) 



TABLE II: Real and imaginary parts of the annihilation contributions for P and S waves. Here /i = rrib/2. 



P-wave 


S'-wave 


/3 2 (7rM 2 ) (0.006; 0.0007) 


(0.031; 0.013) 


/3 3 (ttM 2 ) (-0.024; -0.011) 


(0.094; 0.051) 


k,Ew{-xM 2 ) (0.025; 0.005) x 10~ 2 (- 


-0.009; -0.003) x 10~ 2 



n 

The and B° mean lives, entering in Eq. (|16p . are |5:] t b ± = (1.638 ± 0.011) x 10~ 12 s and 
r B o = (1.530 ± 0.009) x 10~ 12 s, respectively. 

For the Wilson coefficients, C n (//), we take, at both scales \x = m,b and m b /2, the next-to-leading 
order logarithmic approximation values as given in Table 1 of Ref. Qj. Using Eqs. (|18|) and (|19p . 
one obtains, at the scale fi = m b , the universal leading order (LO) ak{^) values presented in the 
first column of Table HI 
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B. Model parameters 

For the B — > ir transition form fact or. \vc> employ tiro pole-extrapolation model [S\. 

if-V) = 7 7^, (55) 

at the momentum transfer, q. In the transition form factor model we are using, the numerical 
parameters are /o(0) = 0.29, a± = 0.76 and 02 = 0.28. 

As pointed out in Sec. HH we use the vector ff K {q 2 ) and scalar ffi K (q 2 ) (with f K /f n = 0.193) 
form factors derived in 

The hard scattering and annihilation contributions for the i^*(892) given in Eqs. ()35[) and (|36|) 
as well as those for the Kq(1430) given in Eqs. (0S]) and (|5"0"|) involve divergences, Xh and Xa 
which are modeled as follows, 

Xa,h = (l + PA,H e-x.p{i(/)A,H)J In (56) 

with, for each X^^, two real parameters pa,h > and < 4>a,h < 360°. One expects the 
annihilation and hard scattering contributions to be of the order of lnfrnfl/Aj,) with = 0.5 GeV 
(see Ref. fl). 

VI. RESULTS AND DISCUSSION 

Within the QCDF approach including final state interactions, before and after hadronization, 
we fit, with the two complex parameters (pa, 4>a) and (pn, <Ph) the mass and helicity angle 
distributions, the P-wave branching ratios and the CP asymmetries provided by the Belle 



and BABAR [13J416J Collaborations. We consider 206 effective mass distribution data, 82 helicity 
distribution points, 6 values of asymmetries for both 7ri^*(892) and 7^^(1430) and 4 branching 
ratios for 7r J fiT*(892). Altogether we have 298 observables with equal weight. Note that in the fit 
we did suppress some points which lie outside the general trend of the data. We have checked that 
these suppressions do not influence the results of the fit. 

We obtain a x 2 /dof = 492.5/(298-4) = 1.68 with the following values p H = 54.43±7.32, 4> H = 
—0.95 ± 0.10 radians for the hard-scattering parameters and pa = 2.51 ± 0.11, <pA = —2.98 ± 0.06 
radians for the annihilation parameters. The corresponding hard-scattering, H n (M2), contributions 
to the short distance amplitudes, ak{p) of Eq. (fT8|) . are listed in Table U together with the leading 
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FIG. 1: Here ABw represent the contributions, to the numerator of the CP asymmetry parameter Aqp 
of Eq. ([TT]). of the different S, P, S + P amplitudes where the scalar and vector form factors have been 
factorized out. The curves denoted by Bg p s+p correspond to the contribution for the S, P, S + P of this 
weak interaction plus perturbative QCD interaction amplitudes to the charged B decays and those denoted 
by Bg P S+P the contributions to the neutral B decays. 



order, vertex V n {M2) and penguin P„(M2) contributions for the P and S-waves. The resulting 
annihilation amplitudes, f3 n (nM2), are displayed in Table ITTl 

The amplitude, a' (/•*), for n = 4 to 10 are always corrections to the a\(fi). For n = 1 to 8, the 
modulus of the LO contribution is larger than the modulus of the vertex term, itself larger than 
that of the penguin. The modulus of the hard-scattering contribution is in between 25% to 60% of 
the modulus of the LO term. The vertex, penguin and hard-scattering contributions can be seen as 
corrections to the leading order amplitude whereas for n = 10, H n (M2) gives the main contribution 
to the very small amplitude a\ Q (fj,). The moduli of the annihilation terms (see Table ILT|) are of the 
order of those of the vertex or penguin for both P- and 5-waves. 

Since the present work and that of Ref. [3( (see their Table VI) use the same leading and 
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FIG. 2: As in Fig. [T]but only for the S + P amplitudes including the scalar and vector form factor contri- 
butions. 

next-to-leading order parameters, the P-wave vertex and penguin contributions to the ak are quite 
similar. For the 5-wave, there are some differences in these corrections for a\ and a|. These arise 
from the introduction of Gegenbauer moments up to order 3 in Ref. [jj. The moduli of the a\ are 
about 20% smaller than those of Ref. [3] (see their Table I). This reduction comes mainly from the 
hard-scattering contributions. 

The K^ttT effective mass distributions for B° — > 7r~7r + are globally well fitted: for B" 

decay, X B cllc/ dof = L03 ' x| A BAR/ dof = °- 71 and for B ° deca y> xLlle/ dof = 10 > X B ABAR/ dof = 

2.96. In the case of the B^ — > tt^tt^K^ 1 effective mass distributions, one has XBeiie/ do ^ = 2-55, 
^babar/ do ^ = 2.65, the data being not very well reproduced, in particular for the charged B 
decays, below 0.9 GeV. The helicity angle distributions are well fitted for both decays with a 
X 2 /dof of the order of 1. 

All the results on branching ratios and asymmetries are summarized in Table IIII1 For the 
if* (892) branching ratios, 90% of the x 2 / do1 comes from the B° and B° BABAR data. These 
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TABLE III: Branching fractions B [see Eq. (flBf ] in units of 10 -6 and direct CP asymmetries Acp in % 
[Eq. ([XT]) ] averaged over charge conjugate reactions. The values of the model, calculated by the integration 
of the itItck distribution over the m^K range from m™^ 1 — 0.82 to m^ x = 0.97 GeV for the P wave and 
from 1.0 to 1.76 for the S wave are compared to the corresponding Belle and BABAR results given in 
the fourth column for B and fifth column for Ac p. Model uncertainties arise from the phenomenological 
parameter errors obtained through the minimization. The third column gives the model values without the 
phenomenological hard scattering and annihilation contributions. 



2?(decay channel) 


Model 


H n [p n ] = 


K?exp f^min rrri max\ 
° \" L ttK ' " L ttK 1 




Refs. 


B(B- - 


-s- tt-K*° - 




5.82 ±0.15 


2.17 


5.35 ±0.59 
5.98 ±0.75 


6.45 ±0.71 
7.20 ±0.90 


[10] 
[13] 


B{B° - 


> TT+K*- - 


► tx+txK) 


4 50 ± 21 


1.65 


4.65 ±0.77 
6.47 ±0.72 


5.60 ±0.93 
11.70 ± 1.30 


[9] 
[16] 


B(B- - 






12.11 ±0.32 


7.80 


25.92 ±2.45 
17.64 ±3.60 


32 ± 3.02 
24.5 ±5.0 


[10] 
[13] 


B(B a -> tt+K*- -» ir+irK) 


11.05 ±0.25 


7.45 


24.95 ±3.25 
12.19 ±3.26 


30.80 ±4.01 
25.40 ±6.80 


M 

[14] 


Acp (decay channel) 


Model 


H n [p n ]{TTM 2 ) = 





A exp 

^CP 


Refs. 


A C p(B- 


->■ TT-R* 


— > 7T~irK) 


0.89 ±0.23 


1.29 




-14.90 ±6.75 
3.2 ±5.4 


[10] 
[13] 


Acp(B° 


TT+K*- 


— > tt + ttK) 


-0.99 ±3.42 


7.99 




-14± 12 


[14] 


A C p(B- 


TT~Kf - 


-> TT~TXK) 


0.27±0.10 


0.27 




7.60 ±4.66 
3.20 ±4.60 


[10] 
[13] 


A CP {B° 


-> n+K*- 


-> it + ttK) 


0.75 ±0.90 


-0.68 




17.0 ±26 


[14] 



are incompatible with the corresponding ones from Belle. The P-wave experimental branching 
ratios for — > 7r ± 7r are well reproduced whereas our predictions for the S-wave branching 
ratios do not fully agree with those provided by Belle but do agree better with the BABAR data. 
As discussed in details in Ref. Q|, the determination of the B — > 7t.K"q(1430) branching ratios is 
problematic as the K* (1430) resonance is wide and the result is quite model dependent. However, 
within the factorization and quasi two-body hypotheses, the use of a scalar form factor, determined 
with precision from theory and experiments other than those of B decays, makes our 7t.K"q(1430) 
branching ratio predictions well founded. 

It is difficult to draw any firm conclusions from the small asymmetries obtained from our global 
fit for both P and S waves since the experimental data have large uncertainties. We found that, 
if we introduced some factor in the x 2 to increase the weight of the CP asymmetries, as done in 
Ref. [3J, we obtain a fit of equivalent quality with, indeed, Aqp values closer to the central values 
of the experimental analyzes, in particular for neutral B decays. 

The plots on effective mass and helicity angle distributions, almost identical to those published 
in [jj, will not be given here. For the 5- wave and for m^K ^ 0.8 GeV, the effective mass distri- 
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butions, mainly for the charged B decays, are smaller than those of Ref. [3] which could indicate 
some stronger suppression of the -K"q(800) contribution. 

In relation with the direct CP violation asymmetries, we will focus on the differential difference 
of effective mass branching ratio distributions for charge conjugate channels. In Fig. [T] we draw 
d(AB)/dm nK with AB = B(B -)• tt-kK) - B(B vfvf^T) [numerator of A C p, see Eq. (TT7])] for 
the charged and neutral decays and calculated from the S, P, and S + P amplitudes, where the 
strong interaction scalar and vector form factors have been factorized out. Figure [2] illustrates 
these distribution differences for the full S + P amplitude including these form factors. The weak 
interaction plus the strong interaction before hadronization produces S + P distribution differences 
(see Fig. [T]) negative for m n K below ~1 GeV, positive and increasing above. Including the final 
state interaction after hadronization the S + P distributions, as seen in Fig. [2] are enhanced in the 
vicinity of the IT* (892) resonance, that of the charged channel is positive while that of the neutral 
is negative. The positive enhancement at the 1430 resonance for the B° decays is larger than that 
of the B~. 

The denominator of Aqp giving similar contribution for charged and neutral channels, the 
above behavior of the S + P distributions allows us to understand the model values (calculated 
by integrating distributions over the m^K range quoted in Table IIIII caption) for Ac p displayed in 
Table IIH| knowing that the P-wave contribution dominates in the vector resonance region and the 
S"-wave in the scalar one. One can see that a strong final state interaction after hadronization can 
increase the CP asymmetry. 

VII. SUMMARY AND OUTLOOK 

In the present study, we analyze the K* resonance effects on the direct CP violation in the B — > 
ttttK decay channels. We calculate the amplitudes for the B° — > it~it + K Q and B^ — > tt^tt^ 



decays in the QCD factorization framework [2j, |6j at leading order in h.QCD/ m b and at the next- 
to-leading order in a s . In order to do so, we approximate these three-body processes as quasi 
two-body B decays into 7ri , T*(892) and ^^(1430) since these final state K* resonances dominate 
the irK effective mass region below 2 GeV. All the contributions, before hadronization, i.e., from 
vertex, penguin, hard-scattering and annihilation corrections as well as those after hadronization, 
i.e., from the K* meson resonance formation and decay described by the strong interaction scalar 
and vector form factors, are included. We complete the calculation performed in Ref. by adding 
explicitly the hard scattering and annihilation contributions which are however subject to large 
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uncertainties arising from the presence of end-point divergences. These divergences are modeled 
with two complex parameters; they are the sole fitted parameters entering in the present calculation. 
Thus, as compared to Ref. [2j|, our model involves only 4 real phenomenological parameters instead 
of 8 while reproducing equally well the present data. These 4 parameters are then determined 
through a fit to the available data on mass and helicity angle distributions, branching ratios and 
CP asymmetries originating from Belle and BABAR Collaboration measurements. The large 
experimental uncertainties in CP asymmetries do not yield strong constraints. Producing higher 
statistics experimental data seems to us mandatory in order to improve constraints on models. 
Furthermore, it should sort out the present discrepancies between the Belle and BABAR analyses. 

At this stage one cannot conclude that the data is or is not compatible with the Standard 
Model. Yet, the possibility of new physics effects, as, for instance, in the minimal super symmetric 
Standard Model approach studied in Ref. \y\ . cannot be excluded. However, the theoretical basis 
of our model being restricted to next-to-leading order corrections, the phenomenological terms can 
simulate next-to-next-to leading order (NNLO) effects. It could also take into account charming 
penguin contributions. In principle, the NNLO corrections to hard scattering are amenable to 
convergent integrals which can be evaluated Q ; This contribution could reduce the phe- 

nomenological part of our model amplitudes. The long distance charming penguin amplitudes 
such as those arising from intermediate D^D^ states could be important since the branching 
fractions for the transition B — > D^D^ are quite large. However, their contributions cannot 
be calculated in a QCD pertubative framework. Both NNLO corrections and charming penguin 
amplitudes should be included before being able to give firm statement as to wether or not it is 
necessary to introduce new physics to understand the data, but, this is outside the scope of the 
present study. 

In conclusion, from this analysis, we point out the important following aspects. 

• It constitutes a robust state of the art QCD factorization calculation at next-to-leading order 
in the strong coupling constant. In this framework, the strong phase can be generated dynamically. 
However, the mechanism suffers from end-point singularities which are not well controlled. It is 
now apparent that the Cabibbo-Kobayashi-Maskawa matrix is the dominant source of CP violation 
in flavour changing processes in B decays. The corrections to this dominant source coming from 
beyond the Standard Model are not expected to be large. In fact, the main remaining uncertainty 
lies in the factorization approximation which provides an explicit picture in the heavy quark limit. 
It takes into account all the leading contributions as well as subleading corrections to the naive 
factorization. The soft collinear effective theory (SCET) has been proposed as a new procedure for 
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factorization [18]. It allows one to formulate a collinear factorization theorem in terms of effective 
operators where new effective degrees of freedom are involved, in order to take into account the 
collinear, soft, and ultrasoft quarks and gluons. Following such steps should improve further our 
knowledge of B physics and, eventually, hint at contributions from physics beyond the Standard 
model. 

• It illustrates explicitly how the strong final state interaction after hadronization can enhance 
CP violation asymmetries. The variation of the differential difference of effective mass branching 
ratio distribution for charge conjugate channels as a function of the irK invariant mass over the 
whole range of the K*(892) and ATq (1430) resonances shows that mixing resonance effects, as those 
seen in Fig. [2J can be observed within a window of 100-200 MeV. With the new Large Hadron 
Collider (LHC) providing energy and accuracy (small energy bin), we believe that by exploring 
such windows the LHCb Collaboration should be in a position to perfom accurate measurements 
of CP violation in B to ttttK decays. 

• It confirms the advantage of using, as a consequence of QCD factorization, a scalar form factor 
to describe the 7tETq(1430) final state. The 1^(1430) resonance is very wide and its nonresonant 
part is difficult to evaluate. Thus, the determination of the B — > 7t1£q(1430) branching fractions 
within, in particular, the isobar model, leads to large uncertainties. As advocated in Ref. [jj, a 
parametrization with this scalar form factor, precisely constructed from unitary coupled channel 
equations using experimental kaon-pion T-matrix elements together with chiral symmetry and 
asymptotic QCD constraints, should be used in experimental Dalitz plot analysis. 
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